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Abstract 

The paper centres on the evaluation of the function n(0) = N(@)/No, 
that is the normalized number of islands as a function of coverage G £ 
[0, 1], given No initial nucleation centres (dots) having any degree of spa- 
tial correlation. A mean field approach has been employed: the islands 
have the same size at any coverage. In particular, as far as the random 
distribution of dots is concerned, the problem has been solved by con- 
sidering the contribution of binary collisions between islands only. With 
regard to correlated dots, we generalize a method previously applied to 
the random case only. In passing, we have made use of the exclusion 
probability reported in [S. Torquato, B. Lu, J. Rubinstein, Phys.Rev.A 
41, 2059 (1990)], for determining the kinetics of surface coverage in the 
case of correlated dots, improving our previous calculation [M. Tomellini, 
M. Fanfoni, M. Volpe Phys. Rev.B 62, 11300, (2000)]. 

1 Introduction 

For describing the evolution of thin film morphologies several physical processes 
have to be considered, among others the atom condensation and evaporation, 
the adatom surface diffusion, the nucleation and the island growth. During the 
last decades significant advances have been done in modeling, taking advantage 
of both analytical and numerical approaches, the inititial stage of film forma- 
tion, namely the nucleation and the growth processes in the low coverage regime 
[Q, |2] • The methods range from the classical mean field rate equations to the 
scaling theories which are found to be suitable tools for describing both exper- 
imental and Kinetic Monte Carlo results In order to go beyond the early 
stage of the growth, it is mandatory to tackle another physical process, that is 



the collision among the growing clusters. In the case of diffusionlcss clusters 
this process is actually linked to the cluster growth only. As a consequence an 
island which is in general a collection of connected clusters stems from one or 
more nucleation events. As regards the film morphology two mechanisms can be 
distinguished: the impingement and the coalescence. In the former no redistri- 
bution of matter among clusters occurs after a collisional event and the clusters 
retain their individuality. Conversely, in the latter mechanism redistribution of 
matter does occur with conservation of both mass and shape. 

The process of coalescence has been modeled, for the first time, by Vincent 
0] who proposed an analytical solution for the island density kinetics on the 
basis of the Poisson distribution. Briscoe and Galvin developed a statistical 
theory for island coalescence which is in good agreement with the behaviour 
of the island density decays obtained by computer simulations 5 . As far as 
the impingement mechanism is concerned, it has been faced in a certain details 
in ref.jHj. Specifically, the exact solution for the normalized number density 
of islands, as a function of surface coverage n(0) = , where N(Q) is the 
number density of islands at coverage and Nq = A(0), has been determined 
in the form of a series, the so called collision series. This computation holds 
in the case of Poissonian simultaneous nucleation and in the entire range of the 
surface coverage (0 < < 1 ). By means of Monte Carlo simulations it was 
also shown that the first three terms well approximate the series |7j. 

Coalescence and impingement mechanisms of 2D and 3D islands have been 
extensively studied, through Monte Carlo simulations, as a function of the nu- 
cleus shape |H]. This work clearly demonstrates that, independently of the col- 
lision mechanism governing the film formation, the n(0) kinetics is, in fact, the 
same, in spite of the fact that in the time domain the two mechanisms give 
rise to completely different kinetics. As a consequence, the behaviour of rt(0) 
is expected to be the same also in the intermediate cases, namely the partial 
coalescence [Sj. The universal behaviour of the n(0) kinetics can be succesfully 
exploited for determining, from experimental data, the nucleation density at 
saturation. To this end it is sufficient to measure island densities in the high 
coverage regime (0 > 15 — 20%), when islands are sufficiently large and atomic 
resolution microscopy is not required. An application of this approach to the 
growth of both diamond on Si substrate and quaterthiophene films on silica 
substrate has been recently presented in 

Thanks to the universal behaviour of the normalized island density, ?i(0) = 
N j^' , also the non simultaneous nucleation case can be tackled. To this purpose 
one resorts to the classical mean field rate equations for both island and adatom 
densities. These equations can be integrated provided the rate coefficient for 
island collisions is evaluated m . As a matter of fact this coefficient is a function 
of n(0) 0. The universal behavior of this function, together with the key 
role it plays in the rate equation approach to film growth, motivated further 
analytical studies. Recently, on the ground of a mean field approximation a semi- 
analytical approach has been developed which leads to island density kinetics 
that is in excellent agreement with the numerical simulation over the whole 
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range of coverage. The method leads to [T2] . 
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where the function W<j(9) = I In j L (1 — 0)^ 1+, ^g?£ has to be evaluated 
numerically. 

The models discussed so far refer to the case in which nuclei are distributed 
at random throughout the whole surface. However, in nucleation processes 
ruled by adatom diffusion there exists a zone, around each nucleus, where the 
nucleation rate is reduced JT]. In other words, the nucleation events do not 
occur at random on the entire surface and this brings about the establishment 
of a degree of spatial correlation among nuclei. 

One of the motivation of the present article is therefore the extension of 
the aforementioned analytical approach to the more realistic case of spatially 
correlated nuclei. Like in the random distribution of nuclei, also in this case 
rate equations can be employed to deal with non simultaneous nucleation. In 
this instance it is worth noting that rate equations are suitable for describing 
Kinetic Monte Carlo results of correlated nucleation ^2]- However, for want of 
theoretical modeling, in the n(8) function has been determined numerically 
through Monte Carlo simulations (MC) |14j . 

The paper is organised as follows. In the first section we propose a novel 
analytical approach based on the statistical mechanics of rigid disks, for the 
evaluation of the n(O) kinetics in the random case. The rate coefficient for 
island collision is also computed, analytically. In the second section the method 
presented in |12j will be employed to tackle the impingement process in simul- 
taneous nucleation of spatially correlated nuclei. 



2 Results and discussion 
2.1 Random distribution of nuclei 

The following analytical approach for computing n(0), is based on a mean field 
approximation, in the sense that at any given coverage (time), all islands have 
the same shape and appropriate size (disks in the case in point). This means 
that we are actually dealing with a sort of coalescence mechanism, yet, due to 
the universal behaviour of the n(0) function, the result of the computation can 
be applied to the impingement case as well. Collisions involving more than two 
islands will be neglected. 

Let us denote by dP the probability that an island be involved in a binary 
collision. Since we are dealing with binary collision the changing rate of the 
island number is just equal to the rate of "dimer" formation; this rate is equal 
to half the rate at which an island undergoes a collision event. Therefore the 
relation holds 
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dN 1 , 

-N = -2 dP > < 2 > 
where N stands for the number density of islands. A simple closure of eqn|21 
is achieved by setting dP — 2irN(2R)d(2R), R being the radius of the islands, 
that is by considering the radial distribution function of the disks to be equal 
to one. For 9 = NnR 2 , to a first approximation one receives ^ff- ~ —2d<d or 

Eg 

n(Q) = e- 2e . (3) 

As will soon be clear this is a poor description of n(Q) over the entire range 
of coverage, it describes the kinetics in the low coverage regime only. In fact, in 
the limit O — ► the kinetics becomes n(Q) ~ 1 — 29 as already derived in [Hj. 

The reason of this inaccuracy relies in the fact that the radial distribution 
function of a system of impenetrable disks is not equal to one except in the limit 
of large particle separation. A more suitable choice of the probability for binary 
collision is indispensable in order to model the island density kinetics over the 
whole range of coverage. To this end the results achieved on the thermodynamic 
system of rigid spheres ^7], ^H]> E3 can be properly exploited. In particular 
we will make use of the results and the notation of ref . EH ■ 

E v (r) is the probability of finding a region of area irr 2 centered at some 
arbitrary point empty of island centers; E p {r) is the probability that given an 
island (its center) at some arbitrary point, the region of area irr 2 encompassing 
the central island is empty of island centers. The respective density probability 
functions are attained by the derivatives of the exclusion probabilities E v (r) and 
Ep(r): H v (r) = — 9E g^ and H p {r) = — dE Q^ . The expression of the exclusion 
probability E v (r), for a system of particles not necessarely correlated through 
hard core potentials, has been derived in refs. E3> OH ■ 

9n the ground of the definition of the exclusion probability one gets dP = 
H p (a)da, where a = 2R is the average value of the island diameter. The rate 
equation for binary collision then becomes 

^ = -\H p {a)da. (4) 

In the thermodynamic limit, "voids" and "islands" exclusion probabilities 
are linked by the relationship E3 

which holds for r > a. Eqn[S] shows that for a system of hard disks E p (cr) = 1. 
Furthermore in the same limit |21| 



H„ 



p (r) = — (a - - ai) E p (r/a). (6) 
a \ o ' 

where a Q = i±||, fll - and E p (x) = e -e[4a (*=-i)+8 Q1 (x-i)] . 

By using eqnEJin the kinetic equation eqn^J we eventually get 
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dN _ 2(l-g6) JQ _ dQ 

N ~ (1-6)2 r c (9)' [n 

where dQ = jirNd(<j 2 ) is the increment of surface fraction covered by islands, 
a = bi — bo and r c (0) the characteristic collision-time function. The definition 
of the time constant as given by eqn[7|is in fact required in order to treat the non 
simultaneous nucleation by menas of rate equations |22 | .|13 ) . The integration 
of eqn[7]can be performed analitically and gives, 

r>„ 2(l-a)e 

n(6) = (l-9) 2a e 1=*~. (8) 

The validity of our approach has been tested by using Monte Carlo (MC) 
simulations of island density decays in him growth ruled by the impingement 
mechanism g). In fig.l eqn|Hl (for a = 0.564 - 0.128 = 0.436 gj]) has been 
compared to the MC simulation. Remarkably, the agreement between the sim- 
ulation and the analytical result is excellent. In any case also the kinetics for 
a = 0.5 — = 0.5 23] is in very good agreement with the simulation (not shown). 

Another analytical but less precise approach, is based on the exclusion prob- 
ability already derived by us in ref.|24j. By only retaining the Heaviside con- 
tribution in the radial distribution function of hard disks and decoupling the 
multiple integral which gives the argument of the exponential of the E v (r) func- 
tion we found 

E v (r) = e - 7rJVr2 [ 1 +^ 7rJVr2 X[o,„]('-)+^JV^ 2 X(< T ^)('')] ] ( 9 ) 

where Xa( x ) ^ s equal to one if x £ A and is equal to zero if x ^ A. It is worth 
noting that this expression has been successfully employed for modeling the 
kinetics of the surface coverage in case of spatially correlated nucleation [see the 
end of the article and [21]]. From eqn|5]on account of eqn|S]we get 



„ , , d r E v (r)\ a d\nE v (r) 
H p (a) = 



E v (a) dr 
that, once inserted in eqnQ] yields the kinetics 



2ttN<t ( 1 + ^Nncr 2 ) , (10) 



n(9) = e -< 2e+3e >. (11) 

As shown in fig.l this kinetics is in a pretty good agreement with the MC 
simulation over the entire range of coverage |25| . In the same hgure it is also 
displayed the behavior of the solution eqn|3] As appears, the modeling based 
on this last approximation is inadequate for describing the kinetics in the whole 
range of surface coverage; in fact it is in agreement with the MC output only at 
very beginning of the kinetics. 

2.2 Spatially correlated nuclei 

The kinetics of island density in the case of Dirac delta (or heterogeneous) 
correlated nucleation can be determined by the same method as that employed 
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in Although the exclusion probability computed in ref. 20 allows, in 

principle, to treat any kind of correlation, here it is introduced in such a way that 
nucleation in a circular region of radius R} lc around each nucleus is prevented. It 
goes without saying that during the growth the radius of the cluster can exceed 
Rhc- Let us introduce the extended surface, e = ttN R 2 , where R stands for the 
nucleus radius and the quantity 6* = ttN R\ c as a measure of the correlation 
degree of the system. The average value of the edge-to-edge distance among 
islands is given in terms of the exclusion probability as follows 



z{R-Q*) = 



OO 

2 J E v {z + R;<d*)dz 

o 

E V (R;B*) 

OO 

2R J E v [R(l + 

o 

!-©(#;©*) 



(12) 



where in the last equation the identity E V (R;Q*) = 1 — 0(i?;0*) has been 
exploited and £ = z/R; clearly R ^ 0. Therefore the ean!12l reads 

with 

1/2 00 

W 1 (e ] e*) = -^ J E v [R(l + Z);Q*]dZ. (14) 
o 

The average diameter of the island when the fraction of covered surface is 
and N is the island density reads 

e ^ 1/2 



D ^= 2 ^im&j) ' (15) 

while the distance between the centers of the islands can be written as 

d(e, e*) = ; ° (i6) 

^(0,6*) 

where C is a constant to be determined. Apparently, 

d(e,e*) = D(&,e*) + z(Q,e*) (17) 

and, since for = 0*/4, N — No, from ean!17l the C constant is easily deter- 
mined as C = 2W 1 (e*/4; 0*) + v 7 ©*. Using eqnsEH [13 [HandE] we end up 
with 
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n(e,e* 



c-e 1 / 2 



(18) 



_Wi(e,e*)_ 

which holds for 8 > 8*/4. As anticipated, we perform the computation of this 
kinetics by resorting to the exclusion probability of ref . [22] which holds for the 
hard core correlation |26j . As far as this expression is concerned, it is given by 



E v (x;y) = [1 - yx 2 } X [0 , 1/2) (x) + (1 - \)e~ ^^^'^^ (x), 

(19) 

where x — (9 e /9*) 1 / 2 (l + £) and y = 8*. Inserting eqn^5] in ean ll4l the 
numerical integration allows one to determine the decay of the islands as shown 
in the fig. 2 together with the MC simulation. The agreement is satisfactory. 

Before concluding this paper, we take the occasion to apply the exclusion 
probability as derived in ref. [2] to the determination of the kinetics of surface 
coverage in the case of Dirac delta spatially correlated nucleation. In our previ- 
ous article [21] we evaluated the kinetics at hands by using eqnOHfor three values 
of 8*, namely 0.2, 0.7 and 1.5. Although those results were quite satisfactory, 
they can be improved by using the following exclusion probability [2] 



22.(9.; e*) = (i-e e )x [0 , T) +(i - ^) e-^ *)[ 4 §^v^+^i)] X[ ^ ^ 

(20) 

where 7(8*) = 1 ( 1 °lef/^)2 4 ■ The results are displayed in fig. 3 together with 
the MC kinetics which have been described in ref. |24| . 
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3 Figure captions 



Fig.l) Kinetics of island impingement as a function of surface coverage in case 
of simultaneous nuclcation of randomly distributed nuclei. The Monte Carlo 
result and the analytical solution (eqn.8) are shown as dots, and full line, re- 
spectively. The solution cqn.ll is reported as dashed line. The behaviour of 
eqn.3 is displayed as open circles. 

Fig. 2) Kinetics of island impingement as a function of surface coverage in 
case of simultaneous nucleation of spatially correlated nuclei. Symbols are the 
Monte Carlo results for a correlated system according to the hard-core model 
for O* = 0.8. The semi analytical result is shown as full line (eqn.18) 

Fig. 3) Kinetics of surface coverage fraction as a function of the extended 
surface, <d e , for simultaneous nucleation of correlated nuclei according to the 
hard-core model. The Monte Carlo simulations and the analytical results are 
shown as symbols and full lines, respectively. Curves a, b and c refers to 0* = 
0.2, 9* = 0.7, 6* = 1.5, respectively. 
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